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Abstract. In this article, we establish optimality of the Bramble-Pasciak-Xu (BPX) 
norm equivalence and optimality of the wavelet modified (or stabilized) hierarchical ba- 
sis (WHB) preconditioner in the setting of local 3D mesh refinement. In the analysis of 
WHB methods, a critical first step is to establish the optimality of BPX norm equivalence 
for the refinement procedures under consideration. While the available optimality results 
for the BPX norm have been constructed primarily in the setting of uniformly refined 
meshes, a notable exception is the local 2D red-green result due to Dahmen and Kunoth. 
The purpose of this article is to extend this original 2D optimality result to the local 3D 
red-green refinement procedure introduced by Bornemann-Erdmann-Kornhuber (BEK), 
and then to use this result to extend the WHB optimality results from the quasiuniform 
setting to local 2D and 3D red-green refinement scenarios. The BPX extension is reduced 
to establishing that locally enriched finite element subspaces allow for the construction 
of a scaled basis which is formally Riesz stable. This construction turns out to rest not 
only on shape regularity of the refined elements, but also critically on a number of ge- 
ometrical properties we establish between neighboring simplices produced by the BEK 
refinement procedure. It is possible to show that the number of degrees of freedom used 
for smoothing is bounded by a constant times the number of degrees of freedom intro- 
duced at that level of refinement, indicating that a practical implementable version of the 
resulting BPX preconditioner for the BEK refinement setting has provably optimal (lin- 
ear) computational complexity per iteration. An interesting implication of the optimality 
of the WHB preconditioner is the a priori //^-stability of the _L2-projection. The exist- 
ing a posteriori approaches in the literature dictate a reconstruction of the mesh if such 
conditions cannot be satisfied. The theoretical framework employed supports arbitrary 
spatial dimension d > 1 and requires no coefficient smoothness assumptions beyond 
those required for well-posedness in H^. 
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1. Introduction 



In this article, we analyze the impact of local mesh refinement on the stability of mul- 
tilevel finite element spaces and on optimality (linear space and time complexity) of 
multilevel preconditioners. Adaptive refinement techniques have become a crucial tool 
for many applications, and access to optimal or near-optimal multilevel preconditioners 
for locally refined mesh situations is of primary concern to computational scientists. The 
preconditioners which can be expected to have somewhat favorable space and time com- 
plexity in such local refinement scenarios are the hierarchical basis (HB) method [9], the 
Bramble-Pasciak-Xu (BPX) preconditioner [16], and the wavelet modified (or stabilized) 
hierarchical basis (WHB) method [35]. While there are optimality results for both the 
BPX and WHB preconditioners in the literature, these are primarily for quasiuniform 
meshes and/or two space dimensions (with some exceptions noted below). In particular, 
there are few hard results in the literature on the optimality of these methods for vari- 
ous realistic local mesh refinement hierarchies, especially in three space dimensions. In 
this article, the first in a series of two articles [2] on local refinement and multilevel pre- 
conditioners, we first assemble optimality results for the BPX norm equivalence in local 
refinement scenarios in three spacial dimensions. Building on the extended BPX results, 
we then develop optimality results for the WHB method in local refinement settings. 
The material forming this series is based on the first author's Ph.D. dissertation [1] and 
comprehensive presentation of this article can be found in [3, 4, 5, 6]. 

Through some topological or geometrical abstraction, if local refinement is extended 
to d spatial dimensions, then the main results are valid for any dimension d > 1 and for 
nonsmooth PDE coefficients p E Loo{0.). Throughout this article, we consider primarily 
the d = 3 case. But, when the abstraction to generic d is clear, we simply state the 
argument by using this generic d. 

The problem class we focus on here is linear second order partial differential equations 
(PDE) of the form: 



Here, / G L2{Q), p,q E L^oi^), p : Q ^ L(3fJ'^, 3?'^), q : Q ^ "R, where p is a 
symmetric positive definite matrix function, and where g is a nonnegative function. Let 
To be a shape regular and quasiuniform initial partition of Vl into a finite number of d 
simplices, and generate 7i, 72, • • • by refining the initial partition using red-green local 
refinement strategies in d = 3 spatial dimensions. Denote as Sj the simplicial linear C° 
finite element space corresponding to Tj equipped with zero boundary values. The set of 
nodal basis functions for Sj is denoted by ^^^^ = {^p''}^^ where Nj = dim Sj is equal 
to the number of interior nodes in Tj, representing the number of degrees of freedom 
in the discrete space. Successively refined finite element spaces will form the following 
nested sequence: 



Let the bilinear form and the functional associated with the weak formulation of (1.1) 
be denoted as 



We consider primarily the following Galerkin formulation: Find u E Sj, such that 



V ■ {p Vm) + q u = f, u = Oon dfl. 



(1.1) 



So c Si c ... c Sj c ... c H^{n). 




f V dx, u,v E Hq (Q). 



a{u,v) = h{v), Wv E Sj. 



(1.2) 
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The finite element approximation in Sj has the form u^^^ = Yld=i'^i^^P ^ where u = 
{ui, . . . .um^Y denotes the coefficients of u^^'^ with respect to The resulting dis- 
cretization operator A^^^ = {a'{4>k \ 4'\^^)}k]=i must be inverted numerically to deter- 
mine the coefficients u from the linear system: 

A^^^u = F^^\ (1.3) 

where F^-'^ = Our task is to solve (1.3) with optimal (linear) complexity 

in both storage and computation, where the finite element spaces Sj are built on locally 
refined meshes. 

Optimality of the BPX norm equivalence with generic local refinement was shown by 
Bramble and Pasciak [14], where the impact of the local smoother and the local pro- 
jection operator on the estimates was carefully analyzed. The two primary results on 
optimality of the BPX norm equivalence in the local refinement settings are due to Dah- 
men and Kunoth [19] and Bomemann and Yserentant [12]. Both works consider only 
two space dimensions, and in particular, the refinement strategies analyzed are restricted 
2D red-green refinement and 2D red refinement, respectively. In this paper, we extend 
the framework developed in [19] to a practical, implementable 3D local red-green refine- 
ment procedure introduced by Bomemann-Erdmann-Komhuber (BEK) [11]. We will 
refer to this as the BEK refinement procedure. 

HB methods [9, 7, 37] are particularly attractive in the local refinement setting be- 
cause (by construction) each iteration has linear (optimal) computational and storage 
complexity. Unfortunately, the resulting preconditioner is not optimal due to condition 
number growth: in two dimensions the growth is slow, and the method is quite effec- 
tive (nearly optimal), but in three dimensions the condition number grows much more 
rapidly with the number of unknowns [26]. To address this instability, one can employ 
L2-orthonormal wavelets in place of the hierarchical basis giving rise to an optimal pre- 
conditioner [23]. However, the complicated nature of traditional wavelet bases, in partic- 
ular the non-local support of the basis functions and problematic treatment of boundary 
conditions, severely limits computational feasibility. WHB methods have been devel- 
oped [34, 35] as an alternative, and they can be interpreted as a wavelet modification 
(or stabilization) of the hierarchical basis. These methods have been shown to optimally 
stabilize the condition number of the systems arising from hierarchical basis methods on 
quasiuniform meshes in both two and three space dimensions, and retain a comparable 
cost per iteration. 

There are two main results and one side result in this article. The main results es- 
tablish the optimality of the BPX norm equivalence and also optimality of the WHB 
preconditoner — as well as optimal computational complexity per iteration — for the re- 
sulting locally refined 3D finite element hierarchy. Both the BPX and WHB precon- 
ditioners under consideration are additive Schwarz preconditioner s. The BPX analysis 
here heavily relies on the techniques of the Dahmen-Kunoth [19] framework and can be 
seen as an extension to three spatial dimensions with the realistic BEK refinement proce- 
dure [11] being the application of interest. The WHB framework relies on the optimahty 
of the BPX norm equivalence. Hence, the WHB results are established after the BPX 
results. 

The side result is the iJ^-stability of L2-projection onto finite element spaces built 
through the BEK local refinement procedure. This question is currently under intensive 
study in the finite element community due to its relationship to multilevel precondition- 
ing. The existing theoretical results, due primarily to Carstensen [18] and Bramble- 
Pasciak-Steinbach [15] involve a posteriori verification of somewhat complicated mesh 
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conditions after local refinement has taken place. If such mesh conditions are not sat- 
isfied, one has to redefine the mesh. However, an interesting consequence of the BPX 
optimality results for locally refined 2D and 3D meshes established here is i/^ -stability 
of L2-projection restricted to the same locally enriched finite element spaces. This re- 
sult appears to be the first a priori i/^-stability result for L2-projection on finite element 
spaces produced by practical and easily implementable 2D and 3D local refinement pro- 
cedures. 

Outline of the paper. In §2, we introduce some basic approximation theory tools used 
in the analysis such as Besov spaces and Bernstein inequalities. The framework for the 
main norm equivalence is also established here. In §3, we list the BEK refinement condi- 
tions. We give several theorems about the generation and size relations of the neighboring 
simplices, thereby establishing local (patchwise) quasiuniformity. This gives rise to an 
L2-stable Riesz basis in §3.1; one can then establish the Bernstein inequality. In §4, we 
explicitly give an upper bound for the nodes introduced in the refinement region. This 
imphes that one application of the BPX preconditioner to a function has linear (optimal) 
computational complexity. In §5, we use the geometrical results from §3 to extend the 
2D Dahmen-Kunoth results to the 3D BEK refinement procedure by establishing the de- 
sired norm equivalence. While it is not possible to establish a Jackson inequality due to 
the nature of local adaptivity, in §6 the remaining inequality in the norm equivalence is 
handled directly using approximation theory tools, as in the original work [19]. In §7, we 
introduce the WHB preconditioner as well as the operator used in its definition. In §8, 
we state the fundamental assumption for establishing basis stability and set up the main 
theoretical results for the WHB framework, namely, optimality of the WHB precondi- 
tioner in the 2D and 3D local red-green refinements. The results in §8 rest completely on 
the BPX results in §5 and on the Bernstein inequalities, the latter of which rest on the ge- 
ometrical results established in §3. The first a priori i/^ -stability result for L2-projection 
on the finite element spaces produced is established in §9. We conclude in §10. 

2. Preliminaries and the main norm equivalence 

The basic restriction on the refinement procedure is that it remains nested. In other 
words, tetrahedra of level j which are not candidates for further refinement will never 
be touched in the future. Let denote the refinement region, namely, the union of the 
supports of basis functions which are introduced at level j. Due to nested refinement 
fij C rij^i. Then the following hierarchy holds: 

fij c ^]J_l c ■■■ c = (2.1) 

In the local refinement setting, in order to maintain optimal computational complexity, 
the smoother is restricted to a local space Sj, typically 

5/ C Sj C S,, (2.2) 

where Sj := {Ij — Ij-i) Sj and Ij : L2(fi) Sj denotes the finite element interpolation 

operator. Degrees of freedom (DOF) corresponding to Sj and Sj will be denoted by 

and Afj respectively where / stands for fine. (2.2) indicates that A/"/ C //j, typically, //j 
consists of fine DOF and their corresponding coarse fathers. 

The BPX preconditioner (also known as parallelized or additive multigrid) is defined 
as follows: 

X«:=X;2^('^-2) 5^ («,#)#. (2.3) 
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Success of the BPX preconditioner in locally refined regimes relies on the fact the BPX 
smoother acts on a local space as in (2.2). As mentioned above, it acts on a slightly 
bigger set than fine DOF (examples of these are given in [13]). Choice of such a set 
is crucial because computational cost per iteration will eventually determine the overall 
computational complexity of the method. Hence in §4, we show that the overall compu- 
tational cost of the smoother is 0{N), meaning that the BPX preconditioner is optimal 
per iteration. We would like to emphasize that one of the the main goals of this paper, 
as in the earlier works of Dahmen-Kunoth [19] and Bomemann-Yserentant [12] in the 
purely two-dimensional case, is to establish the optimality of the BPX norm equivalence: 
J J 
c,J2^'^UQ,-Qj^iMl < Mm<C2j2^'m{Qj-Qj^i)u\\l, (2.4) 

j=0 j=0 

where Qj is the L2-projection. We note that in the uniform refinement setting, it is 
straight-forward to link the BPX norm equivalence to the optimality of the BPX precon- 
ditioner: 

Ci{Xu,u) < \\u\\jji < C2{Xu,u), 

due to the projector relationships between the Qj operators. However, in the local refine- 
ment scenario the precise link between the norm equivalence and the preconditioner is 
more subtle and remains essentially open. 

The rest of this section is dedicated to setting up the framework to establish the main 
norm equivalence (2.4) which will be formalized in Theorem 2.1 at the end of this sec- 
tion. We borrow several tools from approximation theory, including the modulus of 
smoothness, uJk{f,t,^l)p, which is a finer scale of smoothness than differentiability. It 
is a central tool in the analysis here and it naturally gives rise to the notion of Besov 
spaces. For further details and definitions, see [19, 29]. Besov spaces are defined to be 
the collection of functions / G Lp(Q) with a finite Besov norm defined as follows: 

where the seminorm is given by 

|/|B^,,(fi) := \\{2'^uJkif,'2'\^)p}jehNo\\ig, 

with k any fixed integer larger than s. 

Besov spaces become the primary function space setting in the analysis by realizing 
Sobolev spaces as Besov spaces: 

The primary motivation for employing the Besov space stems from the fact that the char- 
acterization of functions which have a given upper bound for the error of approximation 
sometimes calls for a finer scale of smoothness that provided by Sobolev classes func- 
tions. 

The Bernstein inequality is defined as: 

ujk+i{u,t)p<c{mm{l,t2-^}f\\u\\L,, u e S^, j=0,...,J, (2.5) 

where c is independent of u and j. Usually k = degree of the element and in the case 
of linear finite elements k = 1. Here (3 is determined by the global smoothness of the 
approximation space as well as p. For finite elements, (3 = min{l + r + j^, k + 1}. 
Let 9 J be defined as follows. 

e,,j := sup ^^""'^ff/-' , 9 J := max {1, 9^,j : j = 0, . . . , J} . (2.6) 
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Following [19] we have then 

Theorem 2.1. Suppose the Bernstein inequality (2.5) holds for some real number /3 > 1. 
Then, for each < s < min{/3, 2}, there exist constants < ci, C2 < oo independent of 
u & Sj, J = 0, 1, . . ., such that the following norm equivalence holds: 

J J 
^5^22^||(Q,-Q,_i)«llL < \\u\\l,<c,J2^'^\\iQj-Qj->\\l^^^^J- (2-7) 

J j=o j=o 

Proof See [19, Theorem 4.1]. □ 

We would like to elaborate on the difficulities one faces within the local refinement 
framework. In order Bernstein inequality to hold, one needs to establish that the under- 
lying basis is L2-stable Riesz basis as in (3.8). This crucial property heavily depends 
on local quasiuniformity of the mesh. Hence, Bernstein inequality is established in §5 
through local quasiuniformity and L2-stability of the basis in the Riesz sense. 

A Jackson-type inequality cannot hold in a local refinement setting. This poses a 
major difficulty in the analysis because one has to calculate 9j directly. The missing 
crucial piece of the optimal norm equivalence in (2.7), namely, 9j = 0(1) as J — )■ oo, 
will be shown in (6.12) so that (2.4) holds. This required the operator Qj to be bounded 
locally and to fix polynomials of degree 1 as will be shown in §6. 

3. The BEK refinement procedure 

Our interest is to show optimality of the BPX norm equivalence for the local 3D 
red-green refinement introduced by Bomemann-Erdmann-Kornhuber [11]. This 3D red- 
green refinement is practical, easy to implement, and numerical experiments were pre- 
sented in [11]. A similar refinement procedure was analyzed by Bey [10]; in particular, 
the same green closure strategy was used in both papers. While these refinement pro- 
cedures are known to be asymptotically non-degenerate (and thus produce shape regular 
simplices at every level of refinement), shape regularity is insufficient to construct a sta- 
ble Riesz basis for finite element spaces on locally adapted meshes. To construct a stable 
Riesz basis we will need to establish patchwise quasiuniformity as in [19]; as a result, 
d-vertex adjacency relationships that are independent of shape regularity of the elements 
must be established between neighboring tetrahedra as done in [19] for triangles. 

We first list a number of geometric assumptions we make concerning the underlying 
mesh. Let i7 C 3?^ be a polyhedral domain. We assume that the triangulation Tj of 
f2 at level j is a collection of tetrahedra with mutually disjoint interiors which cover 
Q = [j 

reT. ^- "^^^^ ^'^ generate successive refinements %,Ti, . . . which satisfy the 
following conditions: 

Assumption 3.1. Nestedness: Each tetrahedron (son) r E Tj is covered by exactly one 
tetrahedron (father) r' G TJ-i, and any corner ofr is either a corner or an edge midpoint 
ofr'. 

Assumption 3.2. Conformity: The intersection of any two tetrahedra t,t' E Tj is either 
empty, a common vertex, a common edge or a common face. 

Assumption 3.3. Nondegeneracy: The interior angles of all tetrahedra in the refine- 
ment sequence To, 7i, . . . are bounded away from zero. 

A regular (red) refinement subdivides a tetrahedron r into 8 equal volume subtetrahe- 
dra. We connect the edges of each face as in 2D regular refinement. We then cut off four 
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Subdividing 2-simplices (triangles) 




Figure 1 . 

subtetrahedra at the corners which are congruent to r. An octahedron with three paral- 
lelograms remains in the interior. Cutting the octahedron along the two faces of these 
parallelograms, we obtain four more subtetrahedra which are not necessarily congruent 
to r. We choose the diagonal of the parallelogram so that the successive refinements al- 
ways preserve nondegeneracy [1, 10, 27, 38]. A sketch of regular refinement (octasection 
and quadrasection in 3D and 2D, respectively) as well as bisection is given in Figure 1 . 

If a tetrahedron is marked for regular refinement, the resulting triangulation violates 
conformity A. 3. 2. Nonconformity is then remedied by irregular (green) refinement. In 
3D, there are altogether 2^ = 64 possible edge refinements, of which 62 are irregular. 
One must pay extra attention to irregular refinement in the implementation due to the 
large number of possible nonconforming configurations. Bey [10] gives a methodical 
way of handling irregular cases. Using symmetry arguments, the 62 irregular cases can 
be divided into 9 different types. To ensure that the interior angles remain bounded away 
from zero, we enforce the following additional conditions. (Identical assumptions were 
made in [19] for their 2D refinement analogue.) 

Assumption 3.4. Irregular tetrahedra are not refined further. 

Assumption 3.5. Only tetrahedra r E Tj with L{r) = j are refined for the construction 
ofTj+i, where L{t) = min {j : r E Tj} denotes the level ofr. 

One should note that the restrictive character of A. 3. 4 and A. 3. 5 can be eliminated by 
a modification on the sequence of the tetrahedralizations [10]. On the other hand, it is 
straightforward to enforce both assumptions in a typical local refinement algorithm by 
minor modifications of the supporting datastructures for tetrahedral elements (cf. [22]). 
In any event, the proof technique (see (6.8) and (6.9)) requires both assumptions hold. 
The last refinement condition enforced for the possible 62 irregularly refined tetrahedra 
is stated as the following. 

Assumption 3.6. If three or more edges are refined and do not belong to a common face, 
then the tetrahedron is refined regularly. 

We note that the rf-vertex adjacency generation bound for simplices in which are 
adjacent at d vertices is the primary result required in the support of a basis function so 
that (3.6) holds, and depends delicately on the particular details of the local refinement 
procedure rather than on shape regularity of the elements. The generation bound for 
simplices which are adjacent at d—l,d — 2, . . . vertices follows by using the shape regu- 
larity and the generation bound established for d-vertex adjacency. We provide rigorous 
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generation bounds for all the adjacency types mentioned in the lemmas to follow when 
d = 3. The 2D version appeared in [19]; the 3D extension is as described below. 

Lemma 3.7. Let r and r' be two tetrahedra in Tj sharing a common face f. Then 

\L{r)-L{r')\ < 1. (3.1) 

Proof. If L(r) = L{t'), then < 1, there is nothing to show. Without loss of generality, 
assume that L(r) < L{t'). Proof requires a detailed and systematic analysis. To show 
the line of reasoning, we first list the facts used in the proof: 

(1) L{t') < j because by assumption r' G Tj. Then, L(r) < j. 

(2) By assumption r E Tj, meaning that r was never refined from the level it was 
born L(r) to level j. 

(3) Let r" be the father of r'. Then L(r") = L(r') - 1 < j. 

(4) L{t) < L{t') by assumption, implying L(r) < L{t"). 

(5) By (2), r belongs to all the triangulations from L(r) to j, in particular r G Tl{t"), 
whereby (3) L(r") < j. 

/ is the common face of r and t' on level j. By (5) both r, t" E Tl(t")- Then, A. 3. 2 
implies that / must still be the common face of r and t" . Hence, r' must have been 
irregular. 

On the other hand, L(r) < L{t') — 1 = L{t"). Next, we proceed by eliminating the 
possibility that L(r) < L{t"). If so, we repeat the above reasoning, and t" becomes 
irregular, t" is already the father of the irregular r', contradicting A.3.4 for level L{t"). 
Hence L(r) = L{t") = L{t') — 1 concludes the proof. □ 

By A.3.4 and A.3.5, every tetrahedron at any Tj is geometrically similar to some tetra- 
hedron in To or to a tetrahedron arising from an irregular refinement of some tetrahedron 
in To- Then, there exist absolute constants ci, C2 such that 

ci diam(f) 2'^^^'^ < diam(r) < C2 diam(f) 2"^(^\ (3.2) 

where f is the father of r in the initial mesh. The lemma below follows by shape regu- 
larity and (3.1). 

Lemma 3.8. Let t,t' and ^he tetrahedra in Tj sharing a common edge (two 

vertices) and a common vertex, respectively. Then there exist finite numbers V and E 
depending on the shape regularity such that 

\L{r)-L{T')\<V, (3.3) 

1^^(0-^(01 (3.4) 
Consequently, simplices in the support of a basis function are comparable in size as 
indicated in (3.5). This is usually called patchwise quasiuniformity . Furthermore, it 
was shown in [1] that patchwise quasiuniformity (3.5) holds for 3D marked tetrahedron 
bisection by Joe and Liu [24] and for 2D newest vertex bisection by Sewell [30] and 
Mitchell [25]. Due to the restrictive nature of the proof technique (see (6.8) and (6.9)), we 
focus on refinement procedures which obey A.3.4 and A.3.5. However, due to the strong 
geometrical results available for purely bisection-based local refinement procedures, it 
should be possible to establish the main results of this paper for purely bisection-based 
strategies. 

Lemma 3.9. There is a constant depending on the shape regularity of Tj and the quasi- 
uniformity of To, such that 

^Sc.v...'.7^...n.V«. (3.5, 
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Proof. T and r' are either face-adjacent (d vertices), edge-adjacent {d — 1 vertices), or 
vertex-adjacent, and are handled by (3.1), (3.4), (3.3), respectively. 

< c2l^«-^(^')l^^^(by(3.2)) 
diam(r') diam(r') 

< c 2'^^^^i'^'^> 7^°) (by (3.1), (3.4), (3.3) and quasiuniformity of 7^) 

□ 

3.1. L2-stable Riesz basis. Since patchwise quasiuniformity is estabhshed by (3.5), we 
can now take the first step in establishing the norm equivalence in section 5. In other 
words, our motivation is to form a stable basis in the following sense [29]. 

II X] ~ ||{volume^/^(supp ^p^) Ui}^.gA/-Jz2. (3.6) 

The basis stability (3.6) will then guarantee that the Bernstein inequality (2.5) holds. For 
a stable basis, functions with small supports have to be augmented by an appropriate 
scaling so that ||0-''''||l2{!^) remains roughly the same for all basis functions. This is 
reflected in volume (supp (/>■■'■') by defining: 

Lj^i = min{L(r) : r e Tj, Xi G r}. (3.7) 

Then 

volume(supp 0^^) ~ 2"'^^^''. 

We prefer to use an equivalent notion of basis stability; a basis is called L2-stable Riesz 
basis if: 

II M*0?^||L2{f7) ~ \\{ui]x,eMji2^ (3.8) 

"(7") 

where cpi denotes the scaled basis, and the relationship between (3.6) and (3.8) is given 
as follows: 

= 2'^/2L.',: ^. ^ 2-'^/2^^.- Xi e Mj. (3.9) 

Then (3.8) forms the sufficient condition to establish the Bernstein inequality (2.5). This 
crucial property helps us to prove Theorem 8.2. 

Remark 3.10. The analysis is done purely with basis functions, completely indepen- 
dent of the underlying mesh geometry. Furthermore, our construction works for any 
d-dimensional setting with the scaling (3.9). However, it is not clear how to define face- 
adjacency relations for d > 3. If such relations can be defined through some topological 
or geometrical abstraction, then our framework naturally extends to d-dimensional local 
refinement strategies, and hence the optimality of the BPX and WHB preconditioners can 
be guaranteed in d > 1. One such generalization was given by Brandts-Korotov- 
Krizek in [17] and in the references therein. 

4. Local smoothing computational complexity 
In [11], the smoother is chosen to act on the local space 

5, = span [U{</>?}£.,_..iU{<^'' ^ <l^-''&\ ■ 

Other choices for J\fj are also possible; e.g., DOF which intersect the refinement region 
ilj [2, 14]. The only restriction is that Afj C Vlj. For this particular choice, Afj = {i = 
Nj^i + 1, . . . , Nj} : ^^^^ ^ i = l,..., Nj.i}, the following resuk from [11] 
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establishes a bound for the number of nodes used for smoothing (those created in 
by the BEK procedure) so that the BPX preconditioner has provably optimal (linear) 
computational complexity per iteration. 

Lemma 4.1. The total number of nodes used for smoothing satisfies the bound: 

Y,N,<\nj-'^^N,. (4.1) 

Proof. See [11, Lemma 1]. □ 

A similar result for 2D red-green refinement was given by Oswald [29, page 95]. In the 
general case of local smoothing operators which involve smoothing over newly created 
basis functions plus some additional set of local neighboring basis functions, one can 
extend the arguments from [11] and [29] using shape regularity. 



5. Establishing optimality of the BPX norm equivalence 

In this section, we extend the Dahmen-Kunoth framework to three spatial dimensions; 
the extension closely follows the original work. However, the general case for d > 1 
spatial dimensions is not in the literature, and therefore we present it below. 

For linear g, the element mass matrix gives rise to the following useful formula. 

, / X d+1 d+l 

where, i = 1, . . . , c? + 1 and Xi is a vertex of r, = 2, 3. In view of (5.1), we have that 



l|7(i)||2 _ odL,, volume (supp ff^^ 



Since the min in (3.7) is attained, there exists at least one r G supp (j)f^ such that 
L(r) = Lj^i. By (3.2) we have 

2-^.' ^ diam(r)_ 
diam(r) 

Also, 

E 

volume(supp ~ diam'^(rj), Ti G supp . (5.3) 

j=i 

By (3.5), we have 

diam(rj) ~ diam(r). (5.4) 
Combining (5.3) and (5.4), we conclude 

volume(supp (jy'p) ~ E diam'^(r). (5.5) 
Finally then, (5.2) and (5.5) yield 

2'='^i.» volume (supp 0^^) ~ E 

diam (r) 

is a uniformly bounded constant by shape regularity. One can view the size of any 
tetrahedron in To, in particular size of f , as a constant. The reason is the following: 
A. 3. 4 and A. 3. 5 force every tetrahedron at any Tj to be geometrically similar to some 
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tetrahedron in To or to a tetrahedron arising from an irregular refinement of some tetra- 
hedron in To, hence, to some tetrahedron of a fixed finite collection. Combining the two 
arguments above, we have established that 

U?\\L,in)^l, x.eATj. (5.6) 

Let g = Xlx^eAT, ^j^!'''' ^ 'Sj- For any t e Tj we have that 

\\9\\lir)<c \^n4>?\\lin). (5.7) 

where J\fj^r = {xi G Mj : Xi E r}, which is uniformly bounded in r G TJ and j G 
bA^o- By the scaling (3.9), we get equality in the estimate below. The inequality is a 
standard inverse inequality where one bounds g{xi) using formula (5.1) and by handling 
the volume in the formula by (3.2): 

|^.|2 _ 2-'^^-|(7(xi)p < c2-^^-2'^^-||(?||i^(,). (5.8) 

Now, we are ready to establish that our basis is an L2-stable Riesz basis as in (3.8). This 
is achieved by simply summing up over r G Tj in (5.7) and (5.8) and using (5.6). L2 
stability in the Riesz sense allows us to establish the Bernstein inequality (2.5). 

Lemma 5.1. For the scaled basis (3.9), the Bernstein inequality (2.5) holds for (3 = 3/2 

Proof. (5.6) with (5.7) and (5.8) assert that the scaled basis (3.9) is stable in the sense 
of (3.8). Hence, (2.5) holds by [29, Theorem 4]. Note that the proof actually works 
independently of the spatial dimension. □ 



6. Lower bound in the norm equivalence 
The Jackson inequality for Besov spaces is defined as follows: 

inf \\f-g\\L,<cuM,'^~')p, feLpin), (6.1) 

where c is a constant independent of / and J, and a is an integer. In the uniform re- 
finement setting, (6.1) is used to obtain the lower bound in (2.7). However, in the lo- 
cal refinement setting, (6.1) holds only for functions whose singularities are somehow 
well-captured by the mesh geometry. For instance, if a mesh is designed to pick up the 
singularity at x = of y = then on the same mesh we will not be able to recover a 
singularity at x = lofy = — 1). Hence the Jackson inequality (6.1) cannot hold in 
a general setting, i.e. for / G W^. In order to get the lower bound in (2.7), we focus on 
estimating 6j directly, as in [19] for the 2D setting. 

To begin we borrow the quasi-interpolant construction from [19], extending it to the 
three-dimensional setting. Let r G TJ be a tetrahedron with vertices xi, X2, x^, X4. 

Clearly the restrictions of to r are linearly independent over r where Xi G {xi, X2, X3, X4}. 
Then, there exists a unique set of linear polynomials tpl,-ip2,ipl, ipl such that 

(f)l^\x,y,z)ijjj{x,y,z)dxdydz = Ski, Xk,xi G {xi, X2, X3, X4}. (6.2) 

For Xi G Mj and r E Tj, define a function for Xi E t 



M>^'{x,y,z) = < E.^t"- ■'f^-' >■> V ?y) ; ^ ^^y^ 

* ' ' '0, (x, y, z) ^ supp 
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where Ei is the number of tetrahedra in Tj in supp 0^. By (6.2) and (6.3), we obtain 

{mI^\ 0p^) = / mI^\x, y, z)4)i{x, y, z) dxdydz = 6ki, Xk, xi G Mj. (6.4) 

We can now define a quasi-interpolant, in fact a projection onto Sj, such that 

(g,/)(x,y,z) = if,M^'M'\x,y,z). (6.5) 

As remarked earlier, due to (6.3) the slice operator term Qj — Qj-i will vanish outside 
the refined set ilj defined in (2.1). One can easily observe by (5.6) and (6.4) that 

II A^^'^ ||L2(f^) ~ 1, a;^ G AT,, j G biVo. (6.6) 
Letting = \J{t' E Tj : r fl r' 7^ 0}, we can conclude from (5.6) and (6.6) that 

IIQ./IU.M = II E (/'^/'^)'^i'^IU.M < c||/IU.(n,,.). (6.7) 

We define now a subset of the triangulation where the refinement activity stops, mean- 
ing that all tetrahedra inT*, j < m also belong to Tm- 

r* = {teTj: L{t) < J, VLj^r n r' = 0, Vr' G % with L(r') = j}. (6.8) 

Due to the local support of the dual basis functions m/^^ and the fact that Qj is a projec- 
tion, one gets for g E Sj: 

\\9-Qj9\\L,(r) = 0, teT*. (6.9) 

Since Qj is a projection onto linear finite element space, it fixes polynomials of degree 
at most 1 (i.e. ni(3R'^)). Using this fact and (6.7), we arrive: 

ll^-Qj^l|L2(r) < \\g - P\\L2{r) + \\Qj{P - g)\\L2{r) 

< c\\g-PU,in,.^), rET,\T;. (6.10) 

We would like to bound the right hand side of (6.10) in terms of a modulus of smooth- 
ness in order to reach a Jackson-type inequality. Following [19], we utilize a modified 
modulus of smoothness for / G LpiVt) 

^kif,t,n)i = t-^ [ ||A^/r 

They can be shown to be equivalent: 

The equivalence in the one-dimensional setting can be found in [20, Lemma 5.1]. 
For T a simplex in Sft"' and t = diam(r), a Whitney estimate shows that [21, 28, 33] 

inf ||/-P||Lp(r) <CUk+l{f,t,T)p, (6.11) 

Penfc(5Rd) 

where c depends only on the smallest angle of r but not on / and t. The reason why 
Qj works well for tetrahedralization in 3D is the fact that the Whitney estimate (6.11) 
remains valid for any spatial dimension. \ 7J* is the part of the tetrahedralization Tj 
where refinement is active at every level. Then, in view of (3.5) 

diam(fij-,) ~ 2-\ TETj\Tj. 
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Taking the inf over P E ni(3ft^) in (6.10) and using the Whitney estimate (6.11) we 
conclude 

Wg - Qjg\\L2{T) < cu2{g, 2"-', fij-^)2. 

Recalling (6.9) and summing over r E Tj \ T* gives rise to 

\\g - Qjg\\L2{n) < coj2ig,2'\n)2 < 5^2(5', 2"^fi)2, 

where we have switched from the modified modulus of smoothness to the standard one. 
Since Qj is an orthogonal projection, we have the following: 

\\g - Qjg\\ < \\g - Qjg\\- 

Using the above inequality with (2.6) one then has 

= 0(1), oo. (6.12) 

7. The WHB preconditioner 

In local refinement, HB methods enjoy an optimal complexity of 0{Nj — Nj_i) per 
iteration per level (resulting in 0{Nj) overall complexity per iteration) by only using 
degrees of freedom (DOF) corresponding to iSj. However, HB methods suffer from 
suboptimal iteration counts or equivalently suboptimal condition number. The BPX de- 
composition Sj = Sj^i © (Qj — Qj^i)Sj gives rise to basis functions which are not 
locally supported, but they decay rapidly outside a local support region. This allows for 
locally supported approximations, and in addition the WHB methods [34, 35, 36] can be 
viewed as an approximation of the wavelet basis stemming from the BPX decomposi- 
tion [23]. A similar wavelet- like multilevel decomposition approach was taken in [32], 
where the orthogonal decomposition is formed by a discrete L2-equivalent inner prod- 
uct. This approach utilizes the same BPX two-level decomposition [31, 32]. The WHB 
preconditioner is defined as follows: 

/ft.:=X]2^('^-2)5^K#)#, (7.1) 

where ip^-'^ = {Qj — Qj^ijcjy'p . The WHB preconditioner uses the modified basis (where 
as the BPX preconditioner uses the standard nodal basis) where the projection operator 
used is defined as in (7.5). In the WHB setting, these operators are chosen to satisfy the 
following three properties [5]: 

QAs, = /, (7.2) 

QjQk Qmin{j,k}} 

(7.3) 

11(4- - Qj-i)u^'^\\l2 ~ \\u^'^\\l2, u^'^ E (/, - I,-i)S,. (7.4) 

As indicated in (2.2), the WHB smoother acts on only the fine DOF, i.e. , and 
hence is an approximation to fine-fine discretization operator; A^p^ : — )■ where 
'■= {Qj ~ Qj-i)<Sj and / stands for fine. On the other hand, the BPX smoother acts 
on a slightly bigger set than fine DOF, Afj C J\fj typically, union of fine DOF and their 
corresponding coarse fathers. 

The WHB preconditioner introduced in [34, 35] is, in some sense, the best of both 
worlds. While the condition number of the HB preconditioner is stabilized by inserting 
Qj in the definition of Qj, somehow employing the operators Ij — at the same time 
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guarantees optimal computational and storage cost per iteration. The operators which 
will be seen to meet both goals at the same time are: 

Qk = l[l,+Q-{Ij+i-Ij), (7.5) 

j=k 

with Qj = I. The exact L2-projection Qj is replaced by a computationally feasible 
approximation Q| : L2 — )■ Sj. To control the approximation quality of Qj, a small fixed 
tolerance 7 is introduced: 

UQ'^-QMIl, < iWQjuh,, Wu e L2{n). (7.6) 

In the limiting case 7 = 0, Qfc reduces to the exact L2-projection on Sj by (7.2): 

Qk = Qk h+lQk+l ■ ■ ■ Ij-lQj-l I J = QkQk+l ■ ■ ■ Qj-l = Qk- 

Following [34, 35], the properties (7.2), (7.3), and (7.4) can be verified for Qk as follows: 

• Property (7.2): Let u^''^ e Sk- Since [Ij+i - Ij)u^^^ = and Iju^'"'^ = u^'^) for k < j, 
then [Ij + Q]{Ij+i - Ij)]{u^''^) = u^^\ verifying (7.2) for Qk. It also implies 

Ql = Qk. {1.1) 

• Property (7.3): Let k<l, then by (7.7) 

QkQi = [{h + Qlih+i - h)) . . . {h-i + QUih - ii-i)) Qi]Qi = Qk. (7.8) 
Since QkU G Sk and Sk C Si, then by (7.2) we have 

Qi{Qku) = QkU. (7.9) 

Finally, (7.3) then follows from (7.8) and (7.9). 

• Property (7.4): This is an implication of Lemma 7.1. 

For an overview, we list the corresponding slice spaces for the preconditioners of in- 
terest: 

HB: Sj = {I,~I,^i)Sj, 
BPX: 5/ = (g, -g,_i)5,, 

WHB: 5/ = (Qj - Qj-i)S, = {I - Q,^i)(/, - Ij-i)Sj, Q, as in (7.5). 

The WHB smoother only acts on the fine DOF. Then, in the generic multilevel pre- 
conditioner notation, the WHB preconditioner can be written in the following form: 

J 

Bu := Y^Bff\Q, - Qj^i)u. (7.10) 

Bff is chosen to be a spectrally equivalent operator to fine-fine discretization operator 
A^jij. Since the smoother and property (7.4) both rely on a well-conditioned A^jfj, we 
discuss this next. 

7.1. Well-conditioned A^jj. The lemma below is essential to extend the existing results 
for quasiuniform meshes [34, Lemma 6.1] or [35, Lemma 2] to the locally refined ones. 
Sj'^^ = {Ij — Ij^i)Sj denotes the HB slice space. 



Lemma 7.1. Let Tj be constructed by the local refinements under consideration. Let 

Sj = (/ — Qj-i)Sj be the modified hierarchical subspace where Qj-i is any L2- 
bounded operator. Then, there are constants ci and C2 independent of j such that 

ci\W\\\ < U^Wx < C2U^\\x, X = H\L2, (7.11) 
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holds for any = {I ~ Qj_i)(f)f e Sj with 0^ G S^/\ 

Proof. The Cauchy-Schwarz like inequality [8] is central to the proof: There exists 6 E 
(0, 1) independent of the mesh size or level j such that 

(1-(52)(V0^,V0O < (V(0^ + 0O,V(0^ + 0O), V0^G5,_i,0^G5f^. (7.12) 
(1 - 5^) 110-^11^2 < c|0" + 0^|^i (by Poincare inequality and (7.12)). (7.13) 

Combining (7.12) and (7.13): {l-5'^)\\(pm]jr < + Choosing 0^ = -Qj-i<pf, 

we get the lower bound: (1 - 52)||0/||2^^ < HV'-^H^i. 

Let fij denote the support of basis functions corresponding to J\f-. Due to nested 

refinement, triangulation on is quasiuniform. One can analogously introduce a tri- 

angulation hierarchy where all the simplices are exposed to uniform refinement: := 

{r G 7j : L{t) = j} = Tj\^f. Hence, T/ becomes a quasiuniform tetrahedralization 

and the inverse inequality holds for Sj . To derive the upper bound: The right scaling is 
obtained by father-son size relation, and by the inverse inequalities and L2-boundedness 
of Qj-i, one gets 

< Co2'm^f\\l < c,2'^ (l + ||g,-i|U,)' < c2'^\<pf\\l. 

The slice space Sj'^^ is oscillatory. Then there exists c such that ||(/>-'^||i2 < c2"^-^ ||0-^||^i. 
Hence, < c||(/)-^||^i. The case for X = L2 can be established similarly. □ 

Using the above tools, one can establish that A^jj is well-conditioned. Namely, 

ci2'' < < Xin... < C22'\ (7.14) 

where AJj^j^ and Xj^^^^ are the smallest and largest eigenvalues of A^j, and ci are and C2 
both independent of j. For details see [34, Lemma 4.3] or [35, Lemma 3]. 

8. The fundamental assumption and WHB optimality 

As in the BPX splitting, the main ingredient in the WHB splitting is the L 2 -projection. 
Hence, the stability of the BPX splitting is still important in the WHB splitting. The 
lower bound in the BPX norm equivalence is the fundamental assumption for the WHB 
preconditioner. Utilizing a local projection Qj, BPX lower bound was verified earher 
for 3D local red-green (BEK) refinement procedure. The same result easily holds for the 
projection Qj. Dahmen and Kunoth [19] verified BPX lower bound for the 2D red-green 
refinement procedures. 

Before getting to the stability result we remark that the existing perturbation analysis 
of WHB is one of the primary insights in [34, 35]. Although not observed in [34, 35], 
the result does not require substantial modification for locally refined meshes. Let Cj := 
{Qj — Qj)u be the error, then the following holds. 

Lemma 8.1. Let 7 be as in (7.6). There exists an absolute c satisfying: 

J J 
J2'^'me,\\l<c^'J2^'m{Q,-Q,-i)u\\l, yueSj. (8.1) 

j=0 j=0 



Proof. [34, Lemma 5.1] or [35, Lemma 1]. 



□ 
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We arrive now at the primary result, which indicates that the WHB slice norm is opti- 
mal on the class of locally refined meshes under consideration. 

Theorem 8.2. If there exists sufficiently small 70 such that (7.6) is satisfied for 7 e 
[0, 7o), then 

J 

ll^llwHB = $^22^11(0, -Q,-i)w|lL ~ Ml., u G Sj. (8.2) 
i=o 

Proof. Observe that 

[Qj - Qj-i)u = {Qj - Qj)u - {Qj-i - Qj-i)u + {Qj - Qj^i)u (8.3) 
= ej - ej_i + [Qj - Qj^i)u. 

This gives 

i=0 3=Q j=0 

J 

< c(l + 7') J2 IK^J- - Q3-1ML (using (8.1)) 



3=0 



< 4u\\m 



Let us now proceed with the upper bound. The Bernstein inequality (2.5) holds for 
Sj [1, 19] for the local refinement procedures. Hence we are going to utilize an inequality 
involving the Besov norm || ■ ||^i^ which naturally fits our framework when the moduli 
of smoothness is considered in (2.5). The following important inequality holds, provided 
that (2.5) holds [29, page 39]: 



J 

2 

'2,2 



<cJ22'^\u^^^\\l, (8.4) 

3=0 

for any decomposition such that u = '^'j^qU^-'K u^-'^ G Sj, in particular for u^^^ = 
{Qj — Qj-i)u. Then the upper bound holds due to H^iVt) = i?2 2(^)- CH 

Remark 8.3. The following equivalence is used for the upper bound in the proof of 
Theorem 8.2 on uniformly refined meshes [35, Lemma 4]. 

J 

CiMl.< ^ inf J2'^'mu^'m<cMl.. 

Let us emphasize that the left hand side holds in the presence of the Bernstein inequality 
(2.5), and the right hand side holds in the simultaneous presence of Bernstein and Jack- 
son inequalities. However, the Jackson inequality cannot hold under local refinement 
procedures (cf. counter example in section 6). That is why we can utilize only the left 
hand side of the above equivalence as in (8.4). 

Now, we have all the required estimates at our disposal to establish the optimahty 
of WHB preconditioner for 2D/3D red-green refinement procedures for p G Loo{^)- 
We would like to emphasize that our framework supports any spatial dimension d > 1, 
provided that the necessary geometrical abstractions are in place. 
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Theorem 8.4. IfBPX lower bound holds and if there exists sufficiently small 70 such that 
(7.6) is satisfied for 7 G (0, 70), then for B in (7.10): 

{Bu,u) ~ IImII^i- 

Proof. B^pj is spectrally equivalent to . Since is a well-conditioned matrix, using 
(7.14) it is spectrally equivalent to 2^^/. The result follows from Theorem 8.2. □ 

An extension to multiplicative WHB preconditioner is also possible under additional 
assumptions. These results will not be reported here. 

9. iJ^-STABLE L2-PROJECTION 

The involvement of Qj in the multilevel decomposition makes it the most crucial ele- 
ment in the stabilization. We then come to the central question: Which choice of Qj can 
provide an optimal preconditioner? The following theorem sets a guideline for picking 
Qj. It shows that if ^-stability of the Qj is actually a necessary condition for obtaining 
an optimal preconditioner. 

Theorem 9.1. [34, 35]. If Qj induces an optimal preconditioner, namely for u G S,j, 
-O^llia ~ 11^ nil' ^h^^^ ^^^^^^ absolute constant c such that 

\\Qku\\m < c \\u\\Hi, yk < J. 
Proof. Using the multilevel decomposition and (7.3), we get: 

QkU = Yl^j=o(Qj ~ Qj-i)'^- Since Qj induces an optimal preconditioner, there exist two 
absolute constants ai and 02'. 

J 

ffi||n||ii < 5^22^||(Q, -4_i)M||i^ <(T2||^i||ii, MueSj. (9.1) 

i=o 

Using (9.1) for Qfcw: 

k J 

WQkuWl. < -^2^^11(0, -g,_i)n|lL < -X] 2^^11(0, < -Ml.. 

□ 

As a consequence of Theorem 9. 1 we have 

Corollary 9.2. L2-projection restricted to Sj, Qj\sj '■ L2 Sj, is -stable on 2D and 
3D locally refined meshes by red-green refinement procedures. 

Proof. Optimality of the BPX norm equivalence on the above locally refined meshes was 
already established. Application of Theorem 9.1 with Qj proves the result. Alternatively, 
the same result can be obtained through Theorem 9.1 applied to the WHB framework. 
Theorem 8.2 will establish the optimality of the WHB preconditioner for the local re- 
finement procedures. Hence, the operator Qj restricted to Sj is ii^-stable. Since Qj is 
none other than Qj in the limiting case, we can also conclude the ii^-stability of the 
L2-projection. □ 

Our stability result appears to be the first a priori iJ^ -stability for the L2-projection 
on these classes of locally refined meshes, if ^-stability of L 2 -projection is guaranteed 
for the subset Sj of L2{Vl), not for all of 1/2(0). This question is currently undergoing 
intensive study in the finite element and approximation theory community. The exist- 
ing theoretical results, mainly in [15, 18], involve a posteriori verification of somewhat 
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complicated mesh conditions after refinement has taken place. If such mesh conditions 
are not satisfied, one has to redefine the mesh. The mesh conditions mentioned require 
that the simplex sizes do not change drastically between regions of refinement. In this 
context, quasiuniformity in the support of a basis function becomes crucial. This type of 
local quasiuniformity is usually called as patchwise quasiuniformity. Local quasiunifor- 
mity requires neighbor generation relations as in (3.1), neighbor size relations, and shape 
regularity of the mesh. It was shown in [1] that patchwise quasiuniformity holds also for 
3D marked tetrahedron bisection [24] and for 2D newest vertex bisection [25, 30]. These 
are then promising refinement procedures for which if^-stability of the L2-projection 
can be established. 

10. Conclusion 

In this article, we examined the Bramble-Pasciak-Xu (BPX) norm equivalence in the 
setting of local 3D mesh refinement. In particular, we extended the 2D optimality result 
for BPX due to Dahmen and Kunoth to the local 3D red-green refinement procedure in- 
troduced by Bornemann-Erdmann-Kornhuber (BEK). The extension involved establish- 
ing that the locally enriched finite element subspaces produced by the BEK procedure 
allow for the construction of a scaled basis which is formally Riesz stable. This in turn 
rested entirely on establishing a number of geometrical relationships between neighbor- 
ing simplices produced by the local refinement algorithms. We remark again that shape 
regularity of the elements produced by the refinement procedure is insufficient to con- 
struct a stable Riesz basis for finite element spaces on locally adapted meshes. The 
d-vertex adjacency generation bound for simplices in 'R'^ is the primary result required 
to establish patchwise quasiuniformity for stable Riesz basis construction, and this result 
depends delicately on the particular details of the local refinement procedure rather than 
on shape regularity of the elements. We also noted in §3 that these geometrical proper- 
ties have been established in [1] for purely bisection-based refinement procedures that 
have been shown to be asymptotically non-degenerate, and therefore also allow for the 
construction of a stable Riesz basis. 

We also examined the wavelet modified hierarchical basis (WHB) methods of Vas- 
silevski and Wang, and extended their original quasiuniformity-based framework and 
results to local 2D and 3D red-green refinement scenarios. A critical step in the ex- 
tension involved establishing the optimality of the BPX norm equivalence for the local 
refinement procedures under consideration, as established in the first part of this article. 
With the local refinement extension of the WHB analysis framework presented here, we 
established the optimality of the WHB preconditioner on locally refined meshes in both 
2D and 3D under the minimal regularity assumptions required for well-posedness. An 
interesting implication of the optimality of WHB preconditioner was the a priori H^- 
stability of the L 2 -projection. Existing a posteriori approaches in the literature dictate a 
reconstruction of the mesh if such conditions cannot be satisfied. 

The theoretical framework established here supports arbitrary spatial dimension d > 1, 
and therefore allows extension of the optimality results, the if^-stability of L2-projection 
results, and the various supporting results to arbitrary d > 1. We indicated clearly which 
geometrical properties must be re-established to show BPX optimality for spatial dimen- 
sion d > A. All of the results here require no smoothness assumptions on the PDE 
coefficients beyond those required for well-posedness in H^. 

To address the practical computational complexity of implementable versions of the 
BPX and WHB preconditioners, we indicated how the number of degrees of freedom 
used for the smoothing step can be shown to be bounded by a constant times the number 
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of degrees of freedom introduced at that level of refinement. This indicates that practical 
implementable versions of the BPX and WHB preconditioners for the local 3D refine- 
ment setting considered here have provably optimal (linear) computational complexity 
per iteration. A detailed analysis of both the storage and per-iteration computational 
complexity questions arising with BPX and WHB implementations can be found in the 
second article [2] . 
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